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*5 ' Abstract 

We investigate the influence of static disorder and thermal excitations on excitonic energy trans- 
port in the light-harvesting apparatus of photosynthetic systems by solving the Schrodinger equa- 
tion and taking into account the coherent hoppings of excitons, the rates of exciton creation and 
annihilation in antennas and reaction centers, and the coupling to thermally excited phonons. The 
antennas and reaction centers are modeled, respectively, as the sources and drains which provide 
the channels for creation and annihilation of excitons. Phonon modes below a maximum frequency 
are coupled to the excitons that are continuously created in the antennas and depleted in the 
reaction centers, and the phonon population in these modes obeys the Bose-Einstein distribution 
at a given temperature. It is found that the energy transport is not only robust against the static 
disorder and the thermal noise, but it can also be enhanced by increasing the randomness and tem- 
perature in most parameter regimes. Relevance of our work to the highly efficient energy transport 
in photosynthetic systems is discussed. 
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I. INTRODUCTION 



Photosynthesis is an essential means to obtain energy for all life forms on Earth. Although 
the detailed structures of photosynthetic systems are complicated and species-dependent, 
it is now believed that solar photons are absorbed to produce electronic excited states 
(called excitons) in molecular chromophores that are found in light-harvesting (or antenna) 
complexes. Many antenna complexes, together with reaction center complexes which can 
convert the energy stored in the excitons to chemical energy via biochemical reactions, form 
a transportation network for excitons. The photoexcitation process starts with the exciton 
creation in the antenna complexes, which is followed by exciton transfer across the pigment 
network, and ends with exciton trapping by the reaction center complexes. This energy 
transport process is completed on a 10 - 100 picosecond timescale [l|-|6j]. In recent years 
the high efficiency and quantum coherence of the exciton transfer processes amidst a noisy 



environment have atrracted great interest [7 - |20| . The mechanism 



has been extensively studied with various theoretical models 



21 



or such coherent transport 



261 ] . The exciton transfer 



from antennas to reaction centers is often mode 
considers incoherent hoppings between sites [21 



ed by the semiclassical Forster theory which 
22| | . To account for coherence in the energy 



transport processes, a microscopic description is provided by the Redfield theory where the 
exciton dynamics is solved from a master equation in a reduced space of excitons in the 
weak phonon coupling and Born-Markov approximation 23 [.Based on this theory, Silbey 
et al. developed an app roach for the diffusion of excitons 24], |25[ ; Haken and Strobl used 



a stochastic model 
formalism 28 



27]; and Kenkre and Knox established a generalized master equation 



29] to investigate the coherent and incoherent aspects of the excitonic transfer. 
To make the models more realistic, Zhang et al. introduced a modified Redfield equation to 
include static disorder in the exciton system 30] which was later used to simulate the energy 



31] . Silbey and co-workers 



transfer dynamics in light-harvesting complexes of green plants 
proposed a generalized theory which includes coherent transport effects only within donors 
and acceptors while treating interactions between them with the standard Forster model 



32 



34] . A comparison between Forster, Redfield, and other models is given by Yang and 



Fleming in Ref. 35]. With a master equation approach in the Born-Markov and secular 
approximations, Mohseni et al. carried out computer simulations and found that the energy 
transport efficiency can be enhanced by modulating environmental noise 
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In the energy transfer process in photosynthetic systems, excitons are successionally cre- 
ated in the antennas and then trapped in the reaction centers, generating an uninterrupted 
energy flow from the antennas to the reaction centers. This energy flow is significantly af- 
fected not only by the coherence and decoherence effects during the exciton transfers across 
the pigment network, but also by the creation and annihilation of the excitons in the an- 
tennas and the reaction centers, respectively. In fact, if the excitons are viewed as energy 
carriers in the photosynthetic processes, the creation and annihilation of them in the anten- 
nas and the reaction centers play roles of source and drain of carriers, respectively. Thus, 
the widely used source- network- drain models describing continuous current of electron trans- 
port in quantum systems may also be applicable to treat the energy flow in photosynthetic 
systems if the creation and annihilation of excitons can be properly addressed. 

This work is aimed to describe the excitonic energy transport in photosynthetic systems 
with a full quantum-mechanical source-network-drain model. We include on an equal footing 
effects of scattering, disorder, exciton-phonon interactions, thermal excitations in a noisy 
environment, and efficiencies of antennas and reaction centers. The key issues are how to 
model exciton creations and annihilations with a suitable pair of source and drain, and how 
to include the thermal (temperature) effect in such a quantum- mechanical treatment. We use 
incoming (source) channels and outgoing (drain) channels for the creation and annihilation of 
excitons in antennas and reaction centers, and those channels are semi-infinite with the ends 
connected to different sites of a pigment network so that the exciton current from the source 
to the drain can be sustained. Considering that this process is taking place in a thermal 
bath, we include interactions between excitons and bath phonons in the network. With these 
key ingredients captured by a Hermitian Hamiltonian, the wave functions of exciton and 
phonons can be directly solved via the Schrodinger equation without invoking any classical 
or semiclassical approximations. Structural fluctuations in realistic systems are modeled 
with adding disorders to control parameters in the Hamiltonian. Exciton creation and 
annihilation rates can be adjusted by parameter-tuning of incoming and outgoing channels. 
After obtaining the wave functions, the exciton current, which reflects the efficiency of the 
photosynthesis, is calculated via the extended Landauer-Buttiker formula which works in the 



Hilbert space for many-body states 37H39I]. The numerical results show that the creation 



and annihilation rates in antennas and reaction centers and their adaptability have crucial 
effects on the total efficiency and the degree of coherence in the network. It is interesting 
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that the efficiency may be enhanced by introducing randomness and thermal excitations in 
many parameter regimes. At ambient temperature the efficiency may reach its maximum 
or saturation. This can help explain why photosynthesis is so efficient at room temperature 
and in noisy environments. 

The paper is organized as follows. In the next section we present the basic model and 
formalism. In the third section we discuss the results in the absence of interactions with 
phonons and analyze the effect of static randomness. In section IV we include the interaction 
with phonons and investigate the effect of thermal excitations on the energy transport. A 
brief summary is given in section V. 



II. MODEL AND FORMALISM 



We consider the following Hamiltonian for the exciton transport in a photosynthetic 
energy-transfer system: 

H = H C + H s + H d , (1) 

where H c is the Hamiltonian of a pigment network for exciton transport, such as the the light 
harvesting LH1 and LH2 photosynthetic complexes in a membrane of Rsp. photometricum 



40 



41] and the Fenna-Mathews-Olsen (FMO) protein complex in green sulfur bacteria 42], 



connecting the antennas with the reaction centers. H c can be written as 



(2) 



where a\ is the creation operator of exciton on site i, is the site energy for excitons, and Jij 
is the hopping integral for excitons from site j to site i. One has Jjj = J* { as the Hamiltonian 
is Hermitian. Here a site may label a single chromophore or a cluster of chromophores. For 
the former, J it j originates from nearest-neighbor transfer integrals and Forster dipole-dipole 
interactions between well-separated pigments, both of which are strongly dependent on the 



distance between the pigments 



15 



21 



43 



451 ] . For the latter, is the resonance level to 



host an exciton in the ith cluster of chromophores, while Jj j labels the effective transfer 
integral between ith and jth clusters with a distance dependence that differs from that in 
the Forster theory 32J. H s {Hd) is the tight-binding Hamiltonian for the source (drain), 
describing a semi-infinite virtual chain connected to a site in the pigment network as a path 



4 



of the successional creation (annihilation) of excitons in an antenna (reaction center) by 
radiation (biochemical reaction), 

oo 

H s(d)= Yl Ei 5i ( fl t ai ^- i + a l-i a v)+ £ « a L a '>'']' ( 3 ) 

ieS(V) n=l 

where a semi-infinite chain is labeled by the index i of a network site to which it is connected, 
S (T>) is the set of chains belonging to source (drain), n denotes the position on the chain 
counting from the pigment network, and gi is the nearest-neighbor hopping integral in the 
zth virtual chain. Thus the tight-binding Hamiltonian of the zth chain yields a transition 
exciton band of width 4(?j centered at e^, and both gi and vary from chain to chain as 
the transition bands depend on detailed configurations. The widths of transition bands 
in the virtual chains mimic the energy uncertainty of excitons caused by the timescales of 
the photon-exciton conversion events in the antennas and charge separation events in the 



reaction centers. For systems relevant to photosynthesis |46l . 1471 ]. gt varies from 1 meV 
to 10 meV, and e; from 1 eV to 2 eV. As the events occur in a quantized manner, these 
band widths are intrinsic parameters independent of actual light intensity in antennas or 
actual energy flow in reaction canters. In fact, the widths only set upper limits on the 
exciton currents flowing from the sources to the drains, and actual currents may be much 
smaller than those limits as they are also determined by light intensities and adsorption 



cross sections in antennas and by exciton trapping rates in reaction centers |40j. These 
additional physical ingredients, which can be viewed as external conditions that vary from 
time to time, are not included in the Hamiltonian, but their influences on exciton statistics 
in the sources and the drains will be taken into account in the formulation of the exciton 
current. The combined Hamiltonian may be illustrated by Fig. 1, where the networked 
sites corresponding to the antennas (reaction centers) are represented with green (magenta) 
circles, and the creation (annihilation) of excitons is represented with dashed incoming 
(dotted outgoing) arrows attached to corresponding sites. For a specific system such as FMO, 
a site in the pigment network may be connected to an antenna or a reaction center, and it is 
also possible that it is linked only to adjacent sites in the network. Detailed configurations 
in the source-network-drain model can be determined from a realistic photosynthetic system 
of interest. In this work, however, we are concerned with generic features of energy transfer 
in photosynthesis instead of effects of detailed configurations. We will therefore consider 
a network in two dimensions (2D) with closely packed sites with a distribution of control 



FIG. 1: Schematic for a network of exciton transport. The antenna (reaction center) sites 
are represented with green (magenta) circles which are interconnected forming a network, 

and creation (annihilation) of excitons is represented with dashed incoming (dotted 
outgoing) arrows attached to corresponding sites. The couplings between sites, represented 
with brown links in the network, correspond to hopping integrals J^j in Hamiltonian H c , 
while the incoming and outgoing channels expressed by the arrows depict virtual 
semi-infinite chains in Hamiltonians H s and H^. 

parameters. Static disorder erases structural details but the salient features remain after 
configurational averages. The 2D network to be studied here resembles, for example, the 
RC-LH1-LH2 photosynthetic complexes in a membrane of Rsp. photometricum. 

In the photosynthetic process the excitons are constantly flowing into the network via the 
incoming channels. The incident exciton with energy E in an incoming channel i (i G S) 
can be expressed by a plane wave propagating along channel i towards the network 

oo 

^T(E) = ^M~ME)n]4jO), (4) 

n=0 

where ki(E) = arccos[(£' — ei)/2gi\ is the wave vector and E is distributed within the interval 
[e, — 2<7j, e, + 2gi\. The coefficients in the summation are set to unity implying that the rate 
of incident flow is ^gijh if the transition band in the ith chain is occupied by excitons. The 
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incoming wave propagates around in the network and then either exits via the outgoing 
channels or gets reflected to the incoming channels. The outgoing part and reflected part of 
the wave function are also expressed with plane waves in corresponding channels: 

oo 

^(E) = exp[^(£)n]4j0> for j G V, (5) 

n=0 

and 

oo 

^(E) = exp[zM£W4,nl°> for «" e S ' ( 6 ) 

n=0 

where tij and rut are transmission and reflection amplitudes from incident channel % to 
channel j (g X>) and to channel z' (g 5), respectively. The complete form of the wave 
function for exciton with energy E and entering via channel % is 

j€X> i'GS 

Note that the wave function is normalized according to the incident rate. 

The transmission and reflection amplitudes tij(E) and ru/(E) can be solved from the 
Schodinger equation with the transfer-matrix or Green function technique, 

H ^(E) = E^{E). (8) 

The total exciton current through the system, representing the total efficiency of the pho- 
tosynthesis, can be calculated by the Landauer-Buttiker formula summing over all incident 



and output channels and integrating over energy 37j as: 

\tij{E)\ 2 rjij(E)vj(E) \tji(E)\ 2 r]ji(E)vi(E) 



- E 
h ^ 



dE 



Vi(E) Vj(E) 



(9) 



Here \tij(E)\ 2 is the transmission probability from channel i to channel j for given energy 
E, Vi(E) is the exciton velocity in the ith chain, and the ratio Vj(E)/vi(E) accounts for the 
difference in transport rate between the input and output channels. The integration range 
Sij is the region of E where both ki(E) and kj(E) are real. The first and second terms in 
the integrand correspond to the currents from source to drain (positive) and from drain to 
source (negative), respectively. The sign of the total current is controlled by rjij(E), the 
probability that the state of energy E in the iih chain has exciton to emit while the state 
in the jth chain is empty and available to adsorb the exciton. So we can express rjij(E) as 
r)ij(E) = pi(E)qj(E), where Pi(E) (qj(E)) is probability of state in the ith (jth) chain being 



occupied (vacant). As the excitons are created (adsorbed) only in source (drain) chains, 
Pi{E) (qj(E)) is nonzero only when i G 5 (j G D). This guarantees that the second term in 
the above formula is always zero and the current is positive (from source to drain). Pi(E) 
as a function of energy is determined by the light intensity and adsorption cross section in 
antennas, while qj(E) depends on the reopening timescale in reaction centers J40]. Although 
Pi(E) and qj{E) set a further limitation on the energy transport, their effect can be taken 
into account simply by introducing a prefactor to the current. For the sake of simplicity, we 
will adopt a simple form of rjij(E): 

, 1, for i G 5 and 7 G D, 

r)ij{E) = { (10) 
I 0, otherwise. 

Then the current can be calculated as 

I= h P £ dE 1«* sin ' 

where \gj sin kj(E)\/\gi sin ki(E)\ is the ratio of velocities in channels % and j. 

Taking into account the configurational fluctuations in the network, the antennas, and 
the reaction centers, we introduce disorder in and Jy. Since the inter-site coupling in 
the network decreases rapidly as the inter-site distance increases, only couplings between 
nearest neighbors (NN) are kept. The parameter distributions can then be written as 

J l/w a (d), fo r i e 5(D) and e s(d) - io s ( d )/2 < a < e s{d) + w s{d) /2, 
I 0, otherwise, 

j 1/W, for i,j are NN sites and J - < ^ < ^0 + W/2, 

P\Jij) — \ (13) 

I 0, otherwise. 

Here w s (d) and W quantify disorder in for i G 5(D) and in J^- for NN sites, respectively, 
while e s (d) and Jo are their averages. On the other hand, in order to reduce the number of 
parameters, we set gi = g s for i G 5 and ^ = for z G D. 



III. EFFECT OF STATIC DISORDER IN THE ABSENCE OF INTERACTION 
WITH PHONONS 

In this section we investigate general features of energy transfers from the source to the 
drain by taking into account static disorder in the pigment network which reflects configu- 
rational randomness of the photosynthetic system. We carry out numerical calculations for 
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a network of 400 sites closely packed in a 2D plane. At this stage we do not include interac- 
tions with phonons and other thermal excitations, and focus instead on how the efficiency 
of energy transfer is affected by configurations and disorder of the network. Also randomly 
distributed in this network are the source sites connected to the antennas and the drain 
sites connected to the reaction centers. The density of the drain sites is denoted by p. In 
order to reduce the number of parameters in our model, it is assumed that each site in the 
network is connected either to an antenna or a reaction center, so the probabilities of a site 
being source and drain are 1 — p and p, respectively. By simply removing some incoming or 
outgoing channels, our calculation can be easily extended to include cases where some sites 
in the pigment network are connected to neither antennas nor reaction centers. 

Figure 2 displays the network transmission spectrum, i.e., the total probability of energy 
transfer from the incoming photons to the reaction centers as a function of the photon energy, 
for two different amplitudes of drain disorder. The spectrum has a global resonance peak 
of the source-network-drain system. It is apparent that the width of the resonance peak 
is rather small, reflecting the narrow transition bands in the source and drain channels. 
Detuning in the transition bands between the source and the drain rapidly reduces both 
the peak height and width. In the resonant case, e s = e^, the transmission spectrum is 
almost unaffected by the size of the energetic fluctuations in the drain, as can be seen by 
comparing the central peaks in Figs. 2(a) and 2(b). On the contrary, in the detuning case, 
the energetic disorder helps enhance energy transfers by increasing both the height and width 
of the resonance peak. It is interesting to note that even though the transition bands in the 
source and drain channels are in complete detuning, the transmission peak is only partially 
suppressed, but does not vanish. This means that such a pigment network structure with 
multiple antennas and reaction centers is a highly robust system for photosynthesis. 

Now we investigate how the total exciton current, an indicator of the global energy 
transfer efficiency of the system, is affected when the structure and control parameters of 
the pigment network are changed. In Fig. 3, we show the exciton current as a function of the 
reaction center density p for various values of the trapping rate gd which controls the width 
of transition bands in reaction centers. For a given set of g s and gd, there exists a critical 
reaction-center density p c that maximizes the exciton current, i.e., at p c , the excitons created 
in antennas can be most efficiently transported to the reaction centers. Upon increasing the 
transition band width gd, p c is decreased but the exciton current is increased, because for 
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FIG. 2: Total transmission coefficient as a function of the photon energy, (a) 
Wd = 0.002eV, (b) w d = 0.006eV. Other parameters are: e d = 1.5eV, w s = O.OOleV, 
g d = 0.003eV, g s = O.OOleV, J = O.OleV, W = 0.005eV, and p = 0.3. For a given set of 
parameters, there is only a single peak corresponding to the global resonance of the entire 

source-network-drain system. 
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FIG. 3: Total exciton current as a function of drain density p for different values of g d . 
e s = 6d , w d = 0.002eV, w s = O.OOleV, g s = O.OOleV, J = O.OleV, and W = 0.005eV. For 
all combinations of g s and g d the curve of exciton current exhibits a single maximum at p c , 
pointing to an optimal p for energy transfer efficiency 



a larger reaction rate, a given number of excitons can be collected by a smaller number 
of reaction centers. With the given form of rjij(E) in Eq. ffTUl) . the exciton creation and 
trapping rates in the antennas and the reaction centers may be controlled by g s and g d , 
respectively. Thus results obtained here could be used to explain the fact that some purple 
bacteria develop more expansive antenna systems when cultured under low light conditions, 
as a single maximum in the exciton current for a given set of g s and g d indicates that the 
concentration p may be adapted to the environment for optimal exciton transfer efficiency 



481 ] . Our finding is also in qualitative agreement with the work of Fassioli et al. using a 



master equation approach, in which the efficiency can be optimized by changing the number 
ratio of LH2 and LH1 complexes (cf. Fig. 4 of Ref. 40j). Nevertheless, the maxima in Fig. 3 
are not sharp peaks, pointing again to the wide range of structural suitability for the network 
to conduct photosynthetic energy transport. 

To investigate how the energy transfer efficiency is affected by configurational disorder in 
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our model, we plot in Fig. 4 the dependence of the total exciton current on the amplitudes 
of disorder in the pigment network, the sources, and the drains. It is interesting to note 
that the exciton transfer efficiency can be substantially enhanced by increasing the energetic 
disorder in the antennas and the reaction centers. The major efficiency bottlenecks are the 
low creation and annihilation rates in the antennas and reaction centers that lead to a very 
narrow transmission resonance centered at the exciton energies as shown in Fig. 2. Although 
the restriction of the narrow band on the efficiency can be partially lifted by increasing the 
number of channels, there is still a limitation on the maximum current through the drain. 
By introducing the disorder in the exciton energy, the resonant energies spread over a wider 
range, and consequently, the transmission resonance is widened and the exciton current 



increases. In the work of Plenio and co-workers 



49], the efficiency can be enhanced by 



increasing the environment noise in some parameter regimes which correspond to J12 of 



about 100 cm 1 and a full width at hal 



maximum of the spectral density of the noise from 



10 cm 1 to 50 cm 1 (cf. Fig. 4 of Ref. 49j). This is comparable to Jq ~ 0.01 eV and w s ^) 



from 0.0012 eV to 0.0062 eV in Fig. 4. Computer simulations carried out by Mohseni et al. 
also suggest that random noise in the environment might act ually enhance the efficiency of 



the energy transfer in photosynthesis rather than degrade it [36j. On the other hand, the 
transfer rate in the network is usually much larger than the creation and annihilation rates 
in the antennas and in the reaction centers, i.e., Jo ^> g s and Jo ^> ga, so the excitonic 
quantum coherence can be easily preserved during the exciton transfer across the network, 
and the transfer efficiency has a very weak dependence on the network disorder W (i.e., the 
variance in the NN coupling in the network), as shown by the flat I(W) curve in Fig. 4. 
This helps explain why a noisy environment can still preserve quantum coherence and the 
randomness in antennas and reaction centers even enhances the energy transfer efficiency. 



IV. EFFECT OF THERMALLY EXCITED PHONONS 

Previously, the effect of phonons on the energy transfer process was excluded in the 
discussion. However, in realistic photosynthetic systems, thermally excited phonons are 
always present and are believed to play important roles. In this section we investigate the 
effects of phonons and thermal excitations on the energy transport efficiency The phonon- 
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FIG. 4: Dependence of the total exciton current as a function of the degrees of disorder 
Wd (= w s ) and W. For I(w s = Wd), W = 5meV. For I(W), Wd = 2meV, w s = lmeV, and 

e s = €d- Other parameters are: p = 0.3, e s = 1.5eV, ga = 0.003eV, g s = O.OOleV, and 
Jo = O.OleV. The exciton current is enhanced by energetic fluctuations in the antennas and 
the reaction centers, but it is almost unaffected by disorder in the NN hopping integral of 

the pigment network. 



related Hamiltonian can be written as 111, |50|, |51| 



(14) 



k k i 

where bk ip\) is the phonon annihilation (creation) operator for phonon mode k, is the 
corresponding phonon frequency, and labels the exciton-phonon interaction strength at 
site i for phonon mode k. A wave function in the phonon subspace can be written as 



where 



\*P™) = C m,{n k }\<P{n k }), 
{n k } 



i</w>=n4fto> 



'n k \ 



(15) 
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with {nk} being a set of phonon numbers in all modes, and c m> { nfc } is the coefficient in the 
superposition. In the coupled exciton-phonon system, the wave function in Eq. ([7]) can be 
extended to include the phonon Hilbert subspace 

*,,„,,(£) = *;;>,„,,(£) + £E *™lw,«>< B > + EE *a«w.K>^ (16) 

{n'j J^S {n' k } 



where 



and 



*?{n k }( E ) = E eX P[-^-W}(^)^l4n|0) ® I0{n fc} ), (17) 



n=0 



n=0 

(18) 

Here ^fu nh \(E) is the complete wave function for an exciton entering from channel % clothed 
by a phonon cloud |0{ nfe }), and E now labels the total energy of the exciton and companying 
phonons. The first, second, and third terms in Eq. ffTB]) are the incident, transmitted and 
reflected parts of the wave function in different channels and with different phonon states. 
Given the total energy and the phonon energies the wave vector is calculated as ki t { nk y(E) = 
arccos[(.E — — Ylk n k^k) l^9i\- For an exciton-phonon composite with total energy, E 
ti{n h yi {n an d r i {n k }-j {n }( E ) are transmission and reflection amplitudes from channel 
% with phonon state \4>{n k \) to channel j with phonon state |0r„'i), respectively, which can 
be determined from the Schrodinger equation 

{Ho + H ph )V i>{nk} (E) = EV i>{nh} (E) (19) 

with a much-enlarged Hilbert space. By substituting Eqs. (ITT)) and (fl8|) into Eq. (fT4|) . one 
may obtain terms such as 



oo \ / oo 



c i,l 



exp[ik jA n]a] n ®b\\0} + c jA ^2exp[ik jt2 n]a] n ® &J|0> 



,ra=0 / \n=0 



in output channel j, where %,i(2) = arccos[(-E — e,j — %U\(2) )/%9j], an d c j,i( 2 ) * s ^ ne coefficient 
related to the corresponding transmission amplitudes. Obviously these terms cannot be 
factorized if kj t i ^ kj^- Thus, in general, the exciton and phonon states are entangled in an 
output channel even if the two are not in the incoming states. 

As the creation and annihilation of excitons in the antennas and the reaction centers is 
very slow, it is reasonable to assume that when an exciton is created in an antenna, the 
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phonon equilibrium distribution in the network has recovered from the perturbation by the 
movement of the previous exciton. Thus the probability of the phonon state being \4>{ nk }) 
which accompanies an incident exciton is 

P{n k }(T) = —exp[-^2n k hu) k /k B T}, (20) 

k 

where the partition function is defined as 

Z = 22 exp[- n k hcj k / k B T] . 

{n k } k 

Therefore, by including the interactions with phonons, the total exciton current at temper- 
ature T can be calculated as 

where the integration range £ i { nk y.j { n ' } is the span of the total energy E within which both 
k^ {rik} (E) and k j{n ^{E) are real. 

Because of computational constraints, only a restricted number of phonon states are 
included. In the calculation, we adopt a uniform phonon frequency distribution in interval 
[0, Ucl, which is similar to a simplified form of spectral density in Ref. j^]. We adopt ten 
equally spaced modes in this interval, and the phonon states with energies much greater 
than k B T are not included. The interaction strengths \ k j are assumed to be independent 
of the modes, i.e., A^ = A . Calculated temperature dependence of the exciton current is 
plotted in Fig. 5 for various amplitudes of energetic disorder in the sources and the drains, 
and in Fig. 6 for various strengths of exciton-phonon coupling. It is interesting to note that, 
similar to the aforementioned disorder effect, in most cases by increasing the temperature the 
exciton current is enhanced, and eventually reaches saturation at room temperature. This 
means that the thermal excitations may actually favor energy transport in photosynthetic 
systems. The opposite trend, that the current decreases with increasing temperature, only 
occurs in a very narrow range of coupling strength, i.e.j_6 meV < Aq < 7 meV. This result 



is in agreement with that obtained by Mohseni et al. [36[ who has shown that an effective 
interplay between free Hamiltonian evolution and thermal fluctuations in the environment 
leads to a substantial increase in energy transfer efficiency (from about 70% to 99%). Using 
a phonon mode of 180 cm -1 (~ 22meV) and an exciton-phonon coupling strength of ~ 10 
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FIG. 5: Temperature dependence of the total exciton current for various amplitudes of 
energetic disorder in the sources and the drains. Other parameters are: p = 0.3, 
e s = 1.5eV, e d - e s = 2meV, g d = lmeV, g s = 0.5meV, J = O.OleV, W = 5meV, 
u c = 60meV, and Ao = lmeV. For a given value of exciton-phonon coupling strength, e.g., 
Ao = lmeV, the exciton current increases with the increasing temperature for all 

amplitudes of energetic disorder. 



meV, Plenio and coworkers showed that the coupling; to phonons may significantly increase 
energy transfer efficiency (cf. Fig. 7 of Ref. [l4|). Similar values of the exciton-phonon 
coupling strength can also be found in Fig. 6 where an enhancement of the exciton current 
with increasing temperature is shown. The phonon modes which interact with excitons can 
play the role of scatterers for the exciton transport, impeding the energy transmission, but 
they can also provide alternative paths for exciton transfers, therefore promoting energy 
transport. Our results show that at low temperatures the latter effect is dominant in most 
parameter regimes, but at high temperatures the two effects are in balance, yielding a 
temperature- independent exciton current. 
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FIG. 6: Temperature dependence of the total exciton current for various strengths of 
exciton-phonon coupling. Other parameters are: p = 0.3, e s = 1.5eV, — e s = 2meV, 
Wd = w s = lOmeV, gj = lmeV, g s = 0.5meV, Jo = O.OleV, W = 5meV, and u c = 60meV. 
The exciton current is found to increase with the temperature for all exciton-phonon 
coupling strength except the narrow bracket from 6 meV to 7 meV. 



V. CONCLUSIONS 



In this work a source-network-drain exciton Hamiltonian is introduced to investigate the 
energy transport processes in photosynthetic systems. The excitons are electronic energy 
carriers, and the creation and annihilation of the excitons in antennas and reaction cen- 
ters are described by source and drain channels, respectively. The creation and annihilation 
rates in antennas and reaction centers are characterized by the band widths in corresponding 
channels and by the statistical factors of excitons depending on light intensity, adsorption 
cross section in antennas and the reopening timescale in reaction centers. The static dis- 
order in antennas and reaction centers is taken into account by introducing randomness of 
corresponding parameters in the Hamiltonian, while the effects of thermally excited phonons 
are described by exciton-phonon interactions. As the total Hamiltonian is Hermitian and 
the Schrodinger equation is solved for the combined wave function of excitons and phonons, 
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our treatments here are fully quantum mechanical and the effect of coherence is included 
accordingly. The exciton current, an efficiency indicator of the energy transport, can be 
calculated using the Landauer-Buttiker formula. Our results show that the exciton current 
may be enhanced by increasing the static disorder and increasing the temperature in many 
parameter regimes, suggesting that the coherent energy transport in light-harvesting sys- 
tems is robust against most environmental noises, a conclusion that is in agreement with 
previous findings in the literature. Furthermore, the current approach provides a novel, flex- 
ible platform to investigate complicated interplays among various configurational, quantum 
and thermal factors in the energy transport processes of natural and artificial photosynthetic 
systems. 

Previously, many theoretical efforts were based upon solving the phenomenological Lind- 
blad master equations with built-in dissipation. In this work we are interested in the steady 
state in which excitons flow continuously to the reaction centers, and our results are not 
affected by initial conditions and transient behavior. It is also possible to find the steady- 
state density matrix for a single-particle Lindblad equation. We note that the exciton life 
time is in the order of nanoseconds while it takes hundreds of femtoseconds for an exciton to 
be transferred to a reaction center. This huge difference in timescale indicates that exciton 
decay can be neglected when considering excitonic energy transfers from antennas to reac- 
tion centers, therefore lending support to the Landauer-Buttiker scheme adopted here. Our 
approach assumes the single-exciton picture in which the Schrodinger equation is solved, and 
the applicability of this basic assumption to photosynthetic systems is justified by slow injec- 
tion and fast transport of excitons in the pigment network. In a master-equation approach, 
energy transfer efficiency is measured by quantities such as the exciton trapping probability, 
while in the Landauer-Buttiker scheme it is described by the exciton current. Furthermore, 
the total exciton current also provides a global description of the photosynthetic system 
with multi-channel inputs and outputs and interferences among them. 
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